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Abstract — The  Army  Research  Laboratory  (ARL)  has  recently 
developed  the  Recursive  Side-lobe  Minimization  (RSM) 
technique  (patent  pending).  The  technique  is  integrated  with  a 
standard  back-projection  algorithm  to  form  synthetic  aperture 
radar  (SAR)  images  with  significant  reduction  in  side-lobes  and 
noise.  We  have  achieved  significant  improvements  in  noise 
reduction  by  applying  the  RSM  technique  to  our  Ultra- 
wideband  (UWB)  Synchronous  Impulse  Reconstruction  (SIRE) 
forward-looking  radar.  This  paper  presents  the  application  of 
the  RSM  technique  using  data  from  a  side-looking  airborne 
SAR  system  from  SRI  International.  We  describe  the  RSM 
technique,  the  SAR  data  processing,  compare  the  baseline  and 
RSM  SAR  images,  and  quantify  the  image  quality  in  term  of 
signal-to-noise  ratio  (SNR)  and  the  statistical  distribution  of  the 
image  pixels  to  show  that  significant  improvement  achieved 
using  the  RSM  technique. 

1.  Introduction 

Ultra-wideband  (UWB)  synthetic  aperture  radar  (SAR)  systems  have 
been  used  in  many  applications  to  detect  concealed  targets  [1,2]. 
One  of  the  challenges  for  using  SAR  system  in  this  type  of 
application  is  the  detection  of  small  radar  cross  section  (RCS) 
targets  in  the  presence  of  large  man-made  and  natural  clutter  objects. 
Side-lobes  from  large  RCS  objects  spread  throughout  SAR  imagery 
and  create  a  noise  floor.  Signatures  of  small  RCS  targets  might  be 
obscured  or  embedded  in  this  multiplicative  noise  floor  that  is 
proportional  to  the  RCS  of  surrounding  clutter  objects.  In  this  paper, 
we  present  the  Recursive  Side-lobe  Minimization  (RSM)  technique 
(patent  pending)  [3]  designed  to  suppress  the  multiplicative  noise 
floor  in  SAR  imagery.  Although  the  technique  was  originally 
developed  and  tested  using  data  from  the  SIRE  forward-looking 
radar,  we  show  that  is  also  applicable  for  other  SAR  data  sets  with 
different  configurations,  namely  the  SRI  International  airborne  SAR 
radar.  The  paper  first  presents  the  RSM  technique  that  is  integrated 
with  the  back-projection  image  formation  algorithm  to  form  SAR 
imagery.  The  technique  is  then  applied  to  the  SRI  International 
airborne  SAR  radar  data.  Significant  suppression  of  the 
multiplicative  noise  floor  was  achieved  in  the  application  of  the 
RSM  technique  to  this  data  set. 


II.  The  Recurswe  sidelobe  minimization  technique 

In  this  section,  we  present  the  Recursive  Side-lobe 
Minimization  (RSM)  technique  [3],  which  is  designed  to  suppress 
side-lobes  and  noise  from  SAR  imagery.  In  conventional  signal 
processing  techniques  we  usually  try  to  avoid  non-uniform  and 
sparse  data.  From  basic  digital  signal  processing  theory,  it  is  well- 
known  that  if  a  gap  is  introduced  in  the  time-domain  (missing  data 
points);  the  corresponding  frequency  spectrum  would  suffer  from 
the  large  resulting  side-lobes.  This  is  also  true  for  SAR  processing, 
in  which  missing  data  (gaps  in  aperture)  would  generate  sidelobes  in 
imagery  domain.  One  approach  for  improving  SAR  image  in  such 
situations  is  to  estimate  by  interpolating  the  missing  data  samples  in 
the  aperture  data  domain  [4]. 

However,  we  employed  a  different  approach  for  the  RSM  since  we 
already  had  a  complete  data  set.  To  illustrate  the  approach,  we 
generated  a  SAR  simulation  case.  Fig.  1  shows  a  SAR  image  with 
two  simulated  point  targets.  The  targets  have  severe  side-lobes  due 
to  some  sources  (motion  error,  aperture  aliasing)  that  we  introduced 
in  the  SAR  simulation.  Fig.  2a  shows  a  SAR  image  of  the  two  point 
targets  that  is  formed  using  a  randomly  compressed  aperture,  in 
which  approximately  20%  of  the  number  of  aperture  positions  are 
randomly  removed.  As  expected,  the  compressed  aperture  SAR 
image  has  much  larger  side-lobes  than  that  using  of  the  complete 
aperture  due  to  the  gap  introduced  in  the  aperture.  Next  we  formed 
another  SAR  image  (Fig.  2b)  using  another  randomly  compressed 
aperture.  Since  the  locations  of  the  aperture  samples  to  be  removed 
are  randomized,  the  gap  (missing  data)  pattern  would  be  different 
than  the  previous  realization.  We  want  to  note  the  following  two 
important  factors:  1)  Since  the  main  responses  from  the  targets  are 
coherently  integrated  and  normalized  by  the  image  formation 
process,  the  locations  and  the  peak  values  of  the  targets  in  the 
resulting  SAR  images  are  the  same,  and  2)  The  locations  and 
amplitudes  of  peaks  and  valleys  of  the  side-lobes  change  with  each 
formation  of  the  gaps  in  the  compressed  aperture.  The  reason  is  that 
the  only  pixel  locations  in  the  SAR  image  that  are  coherently 
integrated  are  the  target  locations.  Other  pixel  locations  (with  side¬ 
lobes)  are  incoherently  integrated  and  thus  the  values  and  the 
locations  change  with  each  realization  of  compressive  aperture.  Thus 
applying  the  minimizing  operator  to  every  image  pixel  pair  from  the 
two  resulting  compressive  SAR  images  would  yield  the  minimum 
side-lobe  level  in  the  resulting  image  while  the  target  amplitudes 
would  remain  the  same.  Fig.  2c  shows  the  resulting  SAR  image  after 
the  minimum  operator  is  applied  to  every  pair  of  pixels  of  the  two 
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Figure  1 .  SAR  image  of  two 
point  targets  with 
multiplieative  noise. 


This  is  the  main  eoneept  of  the 


b.  The  side-lobe  level  m  Fig.  2e  is 
while  the  amplitudes  of  the  targets 
remain  the  same.  Many 
realizations  of  eompressive 
aperture  are  generated,  eaeh  with  a 
different  gap  pattern.  Due  to  the 
randomness  of  the  loeations  of  the 
peaks  and  valleys  of  the  side- 
lobes,  the  eombined  SAR  image 
using  all  of  these  realizations  has  a 
very  low  side-lobe  level.  Figure  3 
shows  the  results  at  various 
iterations.  In  the  resulting  image 
at  iteration  50,  the  side-lobes  are 
signifieantly  suppressed  while  the 
responses  of  the  two  targets 
remained  virtually  unehanged. 
SM  teehnique. 


Figure  2.  a)  SAR  image  formed  using  a  randomly  eompressed  aperture,  b) 
SAR  image  formed  using  another  randomly  eompressed  aperture,  and  e) 
Resulting  SAR  image  after  the  minimum  operator  is  applied  to  SAR  images 
of  a)  and  b). 


Iteration  1  Iteration  2  Iteration  3 


Figure  3.  Results  at  various  iterations. 


A  step-by-step  deseription  of  the  RMS  algorithm  follows: 

•  Radar  data  and  its  position  information  are  aequired. 

•  The  radar  aperture  is  formed  in  preparation  for  image 
formation.  The  aperture  eonsists  of  K  elements.  Eaeh  element  in 
the  radar  aperture  ineludes  the  radar  reeeiving  position  information 
(x^(A:),y^(A:),z^(A:)),  the  radar  transmitting  information 
{xj{k),yj{k),zj{k)),  and  the  data  reeord  Sk{t)  that  the  radar  measures  at 
this  loeation.  For  side-looking  radar,  the  aperture  is  usually  a  linear 
path  with  data  measured  along  the  path.  For  the  forward-looking 
radar,  a  2-D  radar  aperture  is  generated  [5];  formed  by  the  physieal 
antenna  array  and  the  forward  motion  of  the  radar.  In  general,  the 
radar  aperture  may  take  a  variety  of  shapes. 

•  The  imaging  grid  is  formed  with  N=IJ  pixels  (/  pixels  in 
the  down-range  direetion  and  J  pixels  in  the  eross-range  direetion). 


Eaeh  pixel  in  the  imaging  grid  is  loeated  at  eoordinate 
(xp(/),yp(z),zp(/)),  where  1  <i<N. 

•  A  random  eompressed  aperture  is  generated  using  the 
radar  aperture  with  K  elements  from  step  2  by  seleeting  only  L 
elements  {L<  K)  from  the  original  aperture  for  the  imaging  proeess. 
The  value  for  L  is 

L  =  pK  0<p<l  (1) 

Aeeordingly,  only  a  subset  of  the  aperture  positions  is  used  for 
image  formation.  The  remaining  (K-L)  aperture  positions  are  simply 
disearded  for  this  realization.  The  typieal  number  that  we  use  for 
our  eonfiguration  is  /?  =  0.8  (i.e.,  80%  of  the  aperture  is  employed 
and  20%  of  the  aperture  is  disearded  at  eaeh  iteration).  The  value  of 
p  that  ean  aehieve  best  result  should  be  examined  and  optimized  for 
eaeh  eonfiguration  of  geometry  and  radar  data  set.  The  seleetion  of  L 
aperture  positions  is  eompletely  random  for  eaeh  realization.  If  Ai 
represents  a  veetor  that  eontains  the  indiees  of  aperture  positions  to 
be  ineluded  in  the  image  formation  proeess  for  realization,  then: 

(2) 

where  is  a  random  number,  1  <ain^<K and  ciim  ^orni^m 

•  The  image  is  formed  using  the  eompressed  aperture 
generated  from  previous  step.  The  eompressed  aperture  derived  from 

with  L  elements  is  then  used  to  form  the  realization  of  the  sub¬ 
image  using  the  baek-projeetion  algorithm  [Ref].  This  results  in  the 
realization  of  the  eompressed  aperture  SAR  image. 

I,={P,(i))  \<i<N  (3) 

•  The  envelope  of  the  image  generated  in  previous  step  is 
eomputed.  The  image  generated  in  step  5  ean  also  be  written  as: 

I,={Py))  (4) 

where  Pij  is  the  f  down-range  profile  from  the  r  realization  sub¬ 
image. 

The  eorresponding  quadrature  eomponent  of  this  image  down-range 
profile  is  eomputed  by  applying  the  Hilbert  transform  filter  to  the  in- 
phase  eomponent 

PHy=Hilbert{Py)  (5) 

The  envelope  of  the  down-range  profile  from  the  realization  of 
the  image  may  be  eomputed  as: 

=  (6) 

The  envelope  of  this  image  is  simply 

I,=[PE,j(i))  \<i<I  \<j<J 

•  An  intermediate  resulting  image  is  eomputed.  The 
minimum  operator  is  applied  to  two  images:  1)  the  intermediate 
result  from  previous  iteration  ,  and  2)  the  image  formed  from 
this  iteration.  For  eaeh  image  pixel,  the  values  of  the  two  images  are 
eompared  and  the  minimum  value  is  seleeted 

Im^  =  min(/p  Im;_j )  2<1<M  (8) 

Where  Imi  is  the  intermediate  resulting  image  at  iteration.  Note 

that  equation  (8)  is  defined  for  2  <  /  <  M.  For  the  first  iteration  {1=1), 
Mo  is  initialized  with  a  very  large  values,  so  that  the  intermediate 
resulting  image  is  Im^  =  min(/j,Imo)  = 

•  The  algorithm  returns  to  step  4  to  form  another  randomly 
eompressed  aperture  and  eontinue  with  the  next  iteration  until  the 

iteration  is  finished. 
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III.  APPLICATION  OF  THE  RSM  TECHNIQUE  TO 
THE  SIDE-LOOKING  AIRBORNE  SARDATA 

Although  the  RSM  algorithm  was  originally  developed  for  the 
vehicle-mounted  forward-looking  radar  eonfiguration,  it  is 
applieable  for  other  imaging  seenarios.  In  this  seetion  we  will  show 


Cross-Range 

Figure  4.  ARL  baseline  SAR  image  formed  using  data 
from  the  SRI  airborne  radar. 


that  the  RSM  algorithm  is  also  very  effeetive  to  improve  the  SNR 
using  data  from  the  side-looking  airborne  SAR  data  from  SRI 
International.  The  SRI  radar  is  a  linear  frequeney  modulation  (FM) 
based  UWB  radar  that  operates  in  the  UHF  band.  Our  objeetive  is  to 


a)  SRI  Baseline  Image  b)  ARL  Baseline  Image 


Mean  Target  to  Background  Ratio:  35.6  dB  Mean  Target  to  Background  Ratio:  35.8  dB 


Figure  5.  Baseline  SAR  image  of  ealibration  target  area:  a) 
SRI  baseline  image,  and  b)  ARL  baseline  image. 

proeess  the  SRI  phase  history  data  and  form  a  baseline  SAR  , 
integrate  the  RSM  teehnique  into  the  baseline  algorithm,  and 
evaluate  the  performanee  improvement.  The  proeessing  steps 
inelude  the  range  eompression  and  the  baekprojeetion  image 
formation. 

The  range  eompression  phase  eonsists  of  the  following  steps: 

•  Design  a  matehed  filter 

•  Apply  a  window  funetion  to  the  matehed  filter  for  side- 
lobe  suppression. 

•  Apply  a  Hilbert  Transform  to  the  reeeived  signal  to  form 
an  analytieal  signal. 

•  Down-shift  the  signal  to  baseband. 

•  Apply  the  proper  matehed  filter  to  eompress  the  signal. 

•  Up-shift  the  signal  to  the  transmitted  eenter  frequeney. 

•  Up-sample  the  signal  in  preparation  to  baek-projeetion. 


We  employ  the  standard  baekprojeetion  teehnique  for  SAR  image 
formation.  However,  in  order  to  keep  the  eross-range  resolution 
eonsistent  throughout  a  long  aperture,  we  divide  the  image  area  in 
small  sub-imaging  grids.  For  eaeh  sub-imaging  grid,  we  only  use  a 


Cross-Range 

Figure  6.  SAR  image  after  the  applieation  of  the  RSM 
teehnique  using  data  from  the  SRI  airborne  radar. 


a)  SRI  Baseline  Image  b)  RSM  Image 


Figure  7.  SAR  image  of  a  smaller  area  (400x400  m):  a) 
Baseline  image,  and  b)  RSM  SAR  image. 

segment  of  the  aperture  for  image  formation.  Using  this  approaeh, 
we  ean  aehieve  two  important  features:  ensure  that  the  eross-range 
resolution  and  the  signal-to-noise  (SNR)  level  almost  eonstant 
aeross  a  long  strip  of  SAR  image. 

Fig.  4  shows  the  ARL  baseline  SAR  image  that  eovers  an 
area  of  approximately  2000x2000  m  and  is  generated  using  the 
mosaie  baekprojetion  algorithm  deseribed  above.  The  ARL  baseline 
image  was  formed  using  a  narrow  integration  angle  (9°)  and 
ealibrated  to  mateh  with  the  resolution  and  amplitude  of  the  SAR 
image  provided  by  SRI.  Although  the  9°  integration  angle  is  not  an 
optimal  integration  angle,  it  is  important  that  we  ean  reproduee  the 
baseline  performanee  that  is  eomparable  to  that  of  the  SRI  image. 
From  the  baseline  we  ean  then  measure  the  improvement  aehieved 
by  the  RSM  teehnique.  Fig.  5a  and  5b  show  the  SRI  baseline  and 
ARL  baseline  SAR  images  from  a  ealibration  target  area.  We  define 
the  image  quality  metrie  to  be  the  ratio  of  the  target  amplitude  and 
the  mean  of  the  baekground.  The  target-to-baekground  ratio  is  35.6 
dB  for  the  SRI  baseline  image  and  35.8  dB  for  the  ARL  baseline 
image,  whieh  shows  that  the  ARL  formed  baseline  image  is  well 
matehed  to  the  SRI  formed  image. 
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Fig.  6  shows  the  resulting  SAR  image  after  the  applieation 
of  the  RSM  algorithm.  Compare  Fig.  4  and  Fig.  6,  it  is  obvious  that 
the  baekground  level  in  the  RSM  image  is  signifieantly  lower  than 
that  of  the  SRI  baseline  image.  For  better  visual  eomparison,  we 
display  and  eompare  images  from  several  zoomed-in  areas.  Figure  7 
shows  the  “before”  and  “after”  SAR  images  of  a  smaller  area 
(400x400  m).  It  is  important  to  note  that  the  amplitudes  of  the 
diserete  objeets  in  the  RSM  SAR  image  are  similar  to  that  from  the 
original  SRI  baseline  image  while  the  baekground  level  and  the 
amplitudes  of  the  noise  speekles  are  signifieantly  redueed  in  the 
RSM  SAR  image.  The  enhaneed  SNR  eontrast  aehieved  using  the 
RSM  algorithm  is  obvious  and  signifieant.  Most  of  the  image  pixels 
in  the  RSM  image  are  even  lower  than  the  floor  (-50dB)  of  the 
display  dynamie  range.  Similarly,  the  improvement  is  illustrated  in 
Fig.  8  from  another  image  area. 


a)  SRI  Baseline  Image  b)  RSM  Image 


Figure  8.  a)  SAR  image  of  another  smaller  area  (400x400 
m):  a)  Baseline  image,  and  b)  RSM  SAR  image. 


Mean  Target  to  Background  Ratio:  35.6  dB  Mean  Target  to  Background  Ratio:  41.2  dB 


Figure  9.  SAR  image  of  ealibration  target  area:  a)  SRI 
baseline  image,  and  b)  ARL  RSM  image.  A  5.6  dB  of  improvement 
was  aehieved  using  the  RSM  algorithm. 


To  quantify  the  improvement,  we  eompare  the  targets  and 
baekground  level  from  the  ealibration  target  area  using  the  same 
metrie  that  we  deseribed  earlier  in  this  seetion.  Fig.  9a  and  9b  show 
the  SRI  baseline  SAR  image  and  the  RSM  SAR  image  of  the 
ealibration  target  area.  The  target-to-baekground  ratios  for  the  SRI 
baseline  image  and  RSM  image  is  35.6  dB  and  41.2  dB, 
respeetively.  In  other  words,  a  5.6  dB  of  SNR  improvement  was 
aehieved  by  the  RSM  algorithm. 

Fig.  10  shows  the  improvement  aehieved  by  the  RSM 
algorithm  in  more  details.  The  x-y  plots  show  the  probability 
distributions  of  the  image  pixels  versus  their  amplitudes.  The  blue 
eurve  represents  the  pixel  distribution  of  the  SRI  baseline  SAR 
image  and  the  red  eurve  represents  the  pixel  distribution  of  the  RSM 
SAR  image.  The  pixels  with  high  amplitudes  represent  the  targets 
and  the  image  pixels  with  lower  amplitudes  represent  the 
baekground  noise.  The  plots  show  that  the  RSM  algorithm  retains 
the  amplitudes  from  the  targets  while  signifieantly  suppresses  the 


baekground  noise.  The  figure  also  shows  that  the  RSM  algorithm 
has  moved  the  noise  floor  of  the  original  image  (about  -50  dB)  down 
by  several  dB. 


Figure  10.  The  exeeedanee  distribution  fiinetion  of  the 
image  pixels  versus  their  amplitudes.  The  plots  show  that  the  RSM 
algorithm  retains  the  amplitudes  from  the  targets  while  signifieantly 
suppresses  the  baekground  noise. 

IV.  SUMMARY 

ARL  has  developed  the  Reeursive  Sidelobe 
minimization  (RSM)  algorithm  that  is  integrated  with  the  standard 
baekprojeetion  image  formation  algorithm  to  form  a  SAR  image 
with  signifieant  reduetion  in  sidelobes  and  noise.  The  RSM 
algorithm  has  been  shown  to  be  very  robust  beeause  of  its  simplieity 
(no  training,  no  parameter  tweaking).  Although  the  RSM  was 
initially  developed  for  the  forward-looking  radar  eonfiguration,  we 
have  demonstrated  its  performanee  using  a  different  data  set  (from 
SRI  airborne  radar)  with  a  different  eonfiguration  (side-looking).  For 
this  airborne  side-looking  eonfiguration,  the  RSM  improvement  was 
5.6  dB  over  the  baseline  algorithm.  It  suppresses  many  types  of 
noise  ineluding  the  most  ehallenging  multiplieative  noise  that  is 
originated  from  the  target  signal. 
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